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In this paper the classical Rayleigh-Som- 
merfeld and Kirchhoff boundary-value 
diffraction integrals are solved in closed 
form for circular apertures and slits illu- 
minated by normally incident plane waves. 
The mathematical expressions obtained 
involve no simplifying approximations and 
are free of singularities, except in the 
aperture plane itself. Their use for numeri- 
cal computations was straightforward and 
provided new insight into the nature of dif- 
fraction in the near zone where the 
Fresnel approximation does not apply. The 
Rayleigh-Sommerfeld integrals were 
found to be very similar to each other, so 
that polarization effects appear to be 



negligibly small. On the other hand, they 
differ substantially at sub-wavelength dif- 
ferences from the aperture plane and do not 
correctly describe the diffracted field as 
an analytical continuation of the incident 
geometrical field. 
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1. Introduction 



Diffraction problems in optics typically involve dis- 
tances from the diffracting screen which are large in 
comparison to the wavelength of light. Accordingly the 
Fresnel and Fraunhofer approximations of the principal 
classical diffraction integrals are well documented, but 
so far no workable expressions have been available for 
computations in the near zone. The aim of the present 
paper is to develop mathematical procedures for the 
latter purpose and use them to study the behavior of 
these integrals in the proximity of plane apertures. 

We begin by citing the classical scalar expressions for 
analyzing optical diffraction by an aperture; namely, 
Kirchhoff 's integral equation [1] 



1 r d (e^Q^\ 
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(1) 



and, alternatively, the Rayleigh-Sommerfeld integral 
equations [2, 3, 4] 
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In these equations, U is a. scalar wave function, 5^ is a 
closed surface containing a plane aperture s4 located in 
the xj -plane of a cartesian coordinate system as indi- 
cated in Fig. 1, P = (x, }^, z) is the point of observation 
(z > 0), Q = (^, 17 , 0) is a point on s4, QP is the distance 
between them, n is the aperture normal pointing in the 
direction of the positive z-axis, and k = 27t/\ is the cir- 
cular wavenumber of monochromatic light with wave- 
length A. 
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P = (x,y,z) 




X 



^0 \^o'>yo'^^o) 



Fig. 1. Basic geometry for plane aperture of width 2w. « = aperture normal, 
O = coordinate origin, Pq = point source, P = point of observation. 



2. Background 

The Kirchhoff and Rayleigh-Sommerfeld integral 
equations (1) and (2) are alternative forms of the theo- 
rem of Helmholtz [5], which expresses Huy gens' prin- 
ciple in terms of a scalar wave function U and its normal 
derivatives without assuming specific attributes of this 
function, except that it is continuous and twice differen- 
tiable with continuous derivatives and obeys the homo- 
geneous wave equation, 



AU-^k^U=0, 



(3) 



theorem by itself is insufficient to provide a unique 
solution, it is necessary to impose additional constraints 
on U by prescribing its boundary values on 5^. 

Kirchhoff considered a ''black screen which neither 
reflects nor transmits lights He assumed, plausibly, that 
in this case the incident field vanishes altogether on the 
opaque portion of the screen and is equal to the unper- 
turbed incident field [/geom(Q) inside the aperture. Thus, 

U{Q) = and ^^^ = 0, when Q ^ ^, (4a) 
an 

U(Q) = ^geom(Q), when Q G ^, (4b) 



on and within the closed surface 5^. As Helmholtz' 
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so that Eq. (1) is reduced to Kirchhoff s familiar for- 
mula, where the integration extends over the aperture 
area M only, and U is replaced by f/geom in the integrand. 
As it turned out, Kirchhoff 's solution is mathematically 
flawed. Poincare [6] discovered that it contradicts itself 
and predicted that it will not reproduce the assumed 
boundary conditions Eq. (4a,b). 

Sommerfeld [4, 7] recognized that these difficulties 
are due to the fact that U and b Uldn cannot both vanish 
on any finite portion of the closed surface 5^ unless ^is 
everywhere identically equal to zero. He remedied the 
problem by deriving the integral equations [Eq. (2)] 
which require only the boundary values of either U or 
dU/dn to specify a solution. Thus he assumed, instead 
of Eq. (4a), 



dU(Q) 
dn 



= or U(Q) = 0, when Q ^ 



^(Q) = ^geom(Q), when Q G ^, 



(5a) 



(5b) 



not contribute to the integrals. According to an analysis 
performed by Mukunda [8], these expressions do re- 
cover the assumed boundary conditions in Eqs. (5a-c) as 
P ^ Q, in the sense that U^i replicates the assumed 
value of ^geom but not necessarily the compatible value 
of dUgeom^dn, and the converse is true for U^i. 

It should also be noted that Kirchhoff 's solution is 
simply the arithmetic mean of the Ray leigh- Sommerfeld 
solutions. 



u^(F) = hum) + un(P)l 



(6a) 



and that in the so-called Fresnel limit they are all re- 
duced to one and the same expression. For a point source 
Po = (-^0, Jo, ^o), as shown in Fig. 1, and assuming that the 
distances — zo and z are large in comparison to the aper- 
ture width 2w and the wavelength A, and that 
cos(7T — ^o) and cos 6 are essentially equal to 1, one 
finds [8] 



and applied these conditions separately to the first and 
second Eqs. (2) Except for a passing reference to ''shiny 
screens,'' Sommerfeld did not explicitly associate his 
theory with the diffraction of polarized light but merely 
offered the solution appearing in Eq. (5e), below, as a 
mathematically improved alternative to Kirchhoff 's for- 
mula. On the other hand, Ray leigh [2,3] emphasized that 
the screen must be supposed to be a perfect metallic 
reflector so that the boundary conditions Eq. (5a) per- 
tain top- and s -polarized incident Hght, respectively. He 
emphasized, further, that these conditions apply on the 
dark side of the screen only, while on the lit side it must 
be assumed that 



U^UQ) = or 



a^gom(Q) 

dn 



= 0, when Q ^ , 



(5c) 



According to Maxwell's equations, U^^ and dU^'Vdn 
must be continuous at the screen and hence it follows 
that, taken together, the conditions of Eqs. (5a) and (5c) 
stipulate that for either state of polarization U and 
d Uldn must both zero on the dark side of the screen. 
This is the same as Kirchhoff 's boundary condition of 
Eq. (4a), but without mathematical contradictions. The 
final result for the Ray leigh- Sommerfeld integrals is 
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(5d) 



(5e) 



^k(P) - U^ii^) - ^^s^(P) - U,(F) 

ZZo J 



(6b) 



where, as in Kirchhoff 's theory, the screen itself does 



where /o is the radiant intensity of the source and A(Q) 
is a second-order approximation of the path difference 
(PoQ + QP) - (PqO - OP). It may be estimated that 
this approximation is accurate to 1 % or better when 
z > 20w and z > 20A so that it is usually satisfied for 
narrow apertures and short wavelengths; say, 2w = 0.1 
mm and A = 500 nm, as used for pinhole imagery or 
classroom experiments. In these cases the Rayleigh- 
Sommerfeld and Kirchhoff solutions will hardly be 
needed in their rigorous forms, but on the other hand the 
reliability of the Fresnel approximation is doubtful for 
large apertures and wavelengths. It may not be applica- 
ble in the focal planes of fast lenses or in the case of 
large apertures used in radiometry and photometry, es- 
pecially in the infrared and microwave regions. 

Apart from the above, it appears that the behavior of 
the Kirchhoff and Ray leigh- Sommerfeld integrals in the 
aperture plane has not been documented in the literature 
except in two isolated cases. Wolf and Marchand [10] 
derived a closed expression for Kirchhoff s integral U^ 
inside a circular aperture illuminated by a normally 
incident plane wave, using the Maggi-Rubinowicz trans- 
formation [11] of L^K and a stationary-phase approxima- 
tion. As shown in Fig. 2, this expression gives a fair 
indication of an oscillatory behavior of Uj^ in the aper- 
ture plane but exhibits spurious singularities at the aper- 
ture center and rim. Osterberg and Smith [12] found a 
closed expression for U^i at axial points of observation 
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x/w 



Fig. 2. Approximate behavior of Kirchhoff s integral f/x inside a 
circular aperture of diameter 2w = IDA . Computed from Eq. (A16) of 
Ref. [10]. 



behind a circular aperture and confirmed that it does 
represent a continuous extension of the incident field 
into the half space z > 0. On the other hand, it is evident 
from Fig. 3 that the combined field is not continuously 
differentiate at z = 0, showing that ^^s^ still does not 
fully meet the requirements of Helmholtz' theorem. 
There appear to be no known solutions for U'il, but 
based on Mukunda's work it is clear that 
^^s* = 2Uk— U^l is discontinuous in the aperture plane 
because ^k is discontinuous. 

In addition to these publications, the pertinent litera- 
ture also contains a considerable number of papers that 
were, for the most part, intended to "save" Kirchhoff's 
theory in one way or another. For example, Kottler [13] 
regarded Kirchhoff's integral as the rigorous solution of 
a "saltus" problem in which the wave function U has 
prescribed discontinuities rather than boundary values at 



2 1 



\uf 




kz 



Fig. 3. Behavior of the Rayleigh-Sommerfeld integral Urs for axial 
points near a circular aperture of diameter 2w = lOA . Computed from 
Eq. (7) of Ref. [12]. 



the aperture screen.^ Kottler' s theory involved the 
above-mentioned Maggi-Rubinowicz transformation, 
and subsequently the singularities inherent in the latter 
led to the revival of a belief that diffraction can be 
attributed to "boundary diffraction waves" emerging 
from the edges of aperture screens. Marchand and Wolf 
[14, 15] asserted that the inconsistencies of Kirchhoff's 
integral are only apparent and developed a theory in 
which ^K is expressed in terms of vector potentials that 
have singularities even in free space. Born [16] sug- 
gested the possibility that Kirchhoff's formula may be 
only one in a series of successive approximations, and 
Franz [17] re-derived it by an iterative method in which 
the discontinuities of previous solutions are regarded as 
secondary sources of light. All in all, this curious ex- 
change of conjectures has raised more questions than it 
has answered. Most certainly, it has not addressed the 
concerns of laboratory physicists in search of a "best" 
theory for practical application. 



3. IMathematical Expressions and 
Numerical Results 



3.1 General 

In order to analyze the behavior of the Kirchhoff and 
Rayleigh-Sommerfeld integrals in the proximity of aper- 
tures it is necessary to derive usable expressions for 
computations in the near zone. For this purpose and to 
keep the calculations simple, it will be assumed in the 
following that the incident field is a normally incident 
plane wave so that the geometrical field in the aperture 
is given by 

C/geo„(Q) = VEo, ^^^^^ = ikVEo, (7a) 

where £0 denotes irradiance, and Eqs. (5d) and (5e) can 
be written in normalized form as 



unm = t/&'(P)/V£o 






dQ 
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(7b) 



dz ■ 
(7c) 



^ Saltus is Latin for jump or leap . This author prefers the maxim that 
natura non facit saltus . 
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The simplicity of these relationships is a fortunate 
consequence of having assumed a normally incident 
plane wave. Once the above expression for u^i has been 
evaluated, the solution for u^^ follows by differentiation 
with respect to z, and then Kirchhoffs integral [Eq. (1)] 
is co-determined as the arithmetic mean defined in Eq. 
(6a). There is little doubt that other forms of the incident 
field would have led to considerably more complicated 
expressions without adding to the physical significance 
of the results. In the following, Eqs. (7b) and (7c) will 
be reduced to single integrals for the respective cases of 
circular apertures and slits. 



l3^kQP = kQ^P = k^{^-xf + z^ 



(8a) 



everywhere on these area elements and the integration 
can be carried out over the points Q^ alone. As also 
indicated in Fig. 4, these area elements are in general not 
fully contained in the aperture and must therefore be 
evaluated as 



dQ = 27rd(^-x)(^-x)(l-;^/7r), 



(8b) 



where 2x is the angle subtended by the obstructed arc 
BQ^B' and is given by 



3.2 Circular Aperture 

Let ABCB'A' be the rim of a circular aperture of 
radius 2w which is centered on the origin O of a carte- 
sian coordinate system, as shown in Fig. 4. As the corre- 
sponding diffraction pattern must be rotationally sym- 
metrical about the z-axis, it will be sufficient to 
consider its variation in the xz -plane and the point of 
observation may be chosen as P = {x, 0, z). The inte- 
grals [Eqs. (7b) and (7c)] may then be reduced to single 
integrals by defining the area elements dQ so that they 
are concentric with the projection Qo = (x, 0, 0) of P 
onto the aperture plane and coincide with the circles 
QBQ^B ' shown in the figure, where Q^ = (^, 0, 0) is the 
right-most point at which these circles intersect the x- 
axis. Accordingly, the phases kQP will be constant and 
equal to 



cos x = 



w^-x^-{^-xf 
2x(^-x) 



(8c) 



or ;^ = or IT, as appropriate, when the right-hand side 
of Eq. (8c) exceeds ±1. Consequently, the integrals 
[Eqs. (7b) and (7c)] can be expressed in the form 



u^x^z)-- 



-iejd(^- 



x){^-x){\-xl^)j (8d) 



u^{x,z) = 



ik dz 



-i' 



= k'z\d(^- x)(i - x)(l - x/^){ « - i ^ (8e) 




Fig. 4. Geometrical notation used for circular apertures. 
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the limits of integration being from to w -\- x when 
X < w and from x — w to x-\-w when x > w . 

For the special case of axial points of observation 
(x = 0) the angle x defined by Eq. (8c) is zero, so that 
Eq. (8d) can be solved in closed form. On substitution of 
ip as a new integration variable and subsequent differen- 
tiation with respect to z , one finds 



um.z) = e''' 



W='\/w^-^z\ 



4tO,z) = e^'^- 



ZQ 



W 



(9a) 



(9b) 



the latter being identical to the above-mentioned expres- 
sion derived by Osterberg and Smith [13]. 

For X ^ 0, the numerical integration methods de- 
scribed in Ref. [18] were used to find the real and 
imaginary parts of Eqs. (8d) and (8e) for a small circular 
aperture of diameter 2w = lOA at the distances 
z = O.OIA, A, and lOA. The results obtained are shown 
in Figs. 5a though 5c and will be discussed in Sec. 4. 

3.3 Slit 

Next we consider a diffracting slit of width 2w , cen- 
tered in the xy -plane of a rectangular coordinate system 
as indicated in Fig. 6. The corresponding diffraction 
pattern will consist of straight bands which are parallel 
to the slit jaws, and thus it will again be sufficient to 
compute its variation in the xz -plane. For a given point 
of observation P = (x, 0, z) and for arbitrary aperture 
points Q = (^, 7], 0), Eq. (7b) becomes 



u^^(x,z) = 



W 00 

i/4 



d(kv)- 



.iJ:V(f-.i)2+^2^..2 



w—x 

= ^ j d(^-xM\P), p = kV(^-xf + z\(lO^) 



where Hq^^ = Jo + iYo is the Hankel function of the first 
kind and zero order [19, 20], Jo and Yo are the corre- 
sponding Bessel functions of the first and second kind, 
and j8 is the same as in Eq. (8a). Hence the solution for 
Mrs is obtained at once by substitution of dHo^ 
dz=- khn[^yp into Eq. (10a), leading to 



ukkx,z) = 



ik dz 



/ 




d(^-x)- 



. (10b) 



Fig. 5. Diffraction profiles for a circular aperture of diameter 
2w= lOA at distances z ranging from O.OIA to lOA. The profiles 
shown in each graph are the squared magnitudes of the Rayleigh- 
Sommerfeld solutions U^i and U^i (labeled p and s), and in addition 
the dotted line in Fig. 5c represents the Fresnel approximation Wp 
(labeled F). 



These expressions were evaluated by numerical integra- 
tion, again using the methods of Ref. [18] and assuming 
2w = lOA, z = O.OIA, A, and lOA. The results obtained 
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Fig. 6. Geometrical notation used for slits. 



are shown in Figs. 7a through c. It should be noted that, 
in spite of the singularities of li^\l3) and ll^i\l3)//3 at 
13 = 0, the computations for z = O.OIA presented no 
problems as long as sufficiently small summation ele- 
ments [A(^ — x) = O.Olw] were used. 



4. Discussion 

The mathematical expressions derived in the previous 
Section proved their worth for practical applications, in 
that the computation of the diffraction profiles plotted in 
Figs. 5 and 7 posed no problems. The results obtained 
were everywhere finite, free from singularities, and pro- 
vided new insight into the nature of diffraction in the 
close proximity of apertures. In spite of the obvious 
differences between the profiles pertaining to circular 
apertures slits, their over- all behavior in the near zone is 
similar so that it may be conjectured that the following 
observations are not restricted to these specific aperture 
forms. ^ 



Although Kirchhoff s solution is not explicitly mentioned here, its 
behavior can easily be deduced as it is the arithmetic mean of the 
Rayleigh-Sommerfeld integrals. 



(1) As was to be expected, urs replicates the assumed 
rectangle functions Eqs. (6a,b) in the limit z ^ 0. How- 
ever, it does not constitute an analytical continuation of 
the incident field into the half space z > because other- 
wise du^i/dz, and thus u^i, would also replicate their 
corresponding boundary values. On the other hand, the 
discontinuities of u^i do not manifest themselves in the 
form of sudden jumps, as might be surmised from the 
''saltus'' interpretation of Kirchhoff s theory. Instead, 
they are oscillatory in nature and reminiscent of the 
manner in which the rectangle functions Eqs. (6a,b) 
might be approximated by a Fourier series. 

(2) When z increases, u^i and u^i gradually converge 
to the Fresnel's integral in Eq. (3). For the aperture 
width lOA assumed in the examples the Fresnel limit is 
expected to be reached when z ~ 100 A , and yet Figs. 5c 
and 7c show that the differences between u^i and u^^ are 
already very small at only one tenth this distance. This 
suggests that, according to the Rayleigh-Sommerfeld 
theory, polarization effects are negligibly small even in 
the near zone. 

Although the mathematical expressions derived in 
this paper will be useful for computations in the near 
zone, it remains unclear which of them ought to be used 
in given cases. In a pragmatic sense this may be a mute 
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Fig. 7. Diffraction profiles for a slit of width 2w = IDA at distances z 
ranging from 0.01 A to lOA. The profiles shown in each graph are the 
squared magnitudes of the Rayleigh-Sommerfeld solutions f/S and 
f/Rs (labeled p and s), and in addition the dotted line in Fig. 7c 
represents the Fresnel approximation wp (labeled F). 



question, because u^i and u^l are so similar in most of 
the near zone that either will be an improvement over the 
Fresnel approximation and it does not matter which is 
used. Therefore, u^l could be (and has been) regarded as 
the preferred solution as it is continuous in the aperture 
plane, or Kirchhoff's solution j/k could be regarded as a 
best compromise as it is the arithmetic mean of u^^ and 
t/Rs. In this author's opinion, these are unfounded 
guesses. The fact of the matter is that assessing the 
physical significance of the Rayleigh-Sommerfeld inte- 
grals requires additional considerations that will be the 
subject of a subsequent paper. 
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